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Abstract 

We present a Monte Carlo renormalisation group study of the SU(2) gauge the- 
ory with two Dirac fermions in the adjoint representation. Using the two-lattice 
matching technique we measure the running of the coupling and the anomalous mass 
dimension. We find slow running of the coupling, compatible with an infrared fixed 
point. Assuming this running is negligible we find a vanishing anomalous dimension, 
7 = —0.03(13), however taking this source of systematic error into account gives a 
much larger range of allowed values, —0.6 < 7 < 0.6. We also attempt to measure 
the anomalous mass dimension using the stability matrix method. We discuss the 
systematic errors aflFecting the current analysis and possible improvements. 



1 Introduction 



Technicolor theories with fermions in higher representations of the gauge group can poten- 
tially provide a dynamical electroweak symmetry breaking mechanism without conflicting 
with electroweak precision data. Minimal Walking Technicolor is an example of such a 
theory, a SU(2) gauge theory with two Dirac fermions in the adjoint representation [Il|2]. 
It is expected from perturbation theory to be in or near to the conformal window, al- 
though any new infrared fixed point (IRFP) is thought to occur at strong coupling and so 
non-perturbative results are necessary. Initial lattice simulations showed some evidence 
of walking dynamics and mapped out the phase diagram of the theory [3l H] . Subsequent 
Schrodinger Functional lattice simulations have indeed found that the gauge coupling runs 
very slowly [SI El El E], more slowly than the perturbative prediction, and a recent sim- 
ulation using Creutz ratios found evidence for backwards running at strong coupling [9], 
although distinguishing between conformal and near-conformal behaviour is an inherently 
difficult task. Moreover, in order to be phenomenologically viable the theory must have a 
large anomalous mass dimension (7 ~ 1) [IHl [HI [12] , and recent work suggests that 7 > 1 
is required [T3] . 

A conjectured all-order beta function [H] predicts 7 = 11/24 ^ 0.458 for this modelj^ 
This value is also consistent with perturbative results in the MS-scheme up to four 
loops [16]. The anomalous mass dimension has been measured non-perturbatively in 
recent lattice studies [51 [HI [H], which have all found lower values. 

In this work we measure the anomalous mass dimension using the Monte Carlo Renor- 
malisation Group (MCRG) two-lattice matching method, a technique which has recently 
been used to investigate theories with many flavours of fermions in the fundamental rep- 
resentation [iHl [19] . We also investigate using the stability matrix MCRG method [2^ , 
which in principle allows the determination of all the critical exponents of a system. We 
consider the evolution of all possible couplings of the system under Wilson Renormalisa- 
tion Group (RG) block transformations, where with each blocking step ultraviolet (UV) 
fluctuations are integrated out. Fixed points are characterised by the number of relevant 
couplings, which have positive scaling dimensions and flow away from the flxed point. 
Irrelevant couplings have negative scaling dimensions, and flow towards the flxed point, 
so that their infrared (IR) values are independent of their UV values. 

2 Two-Lattice Matching Method 

With each RG step, changing the scale by a factor s, irrelevant couplings will flow towards 
the flxed point (FP), and relevant couplings will flow away from it. After a few steps the 

^This prediction supersedes the original all-order conjecture 15^ of 7 = 3/4 for this model. 
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irrelevant couplings should die out, leaving the flow following the unique renormalised 
trajectory (RT). If we can identify two sets of couplings which end up at the same point 
along the RT after the same number of steps, then they must have the same lattice 
correlation lengths, ^ = Since the physical correlation length ^ = should not be 
changed by the RG transform, this means that they both must have the same lattice 
spacing a, or inverse cutoff ~ a. If they end up at the same point, but one takes an 
extra step, then their lattice correlation lengths must differ by a factor s, and hence so 
must their UV cutoffs. 

To identify such a pair of couplings, we need to show that after n and {n — 1) RG steps 
respectively their actions are identical. Explicitly calculating the actions would be com- 
plicated, but instead the gauge configurations themselves can be RG block transformed. 
Showing that the expectation values of all observables on these gauge configurations agree 
is equivalent to directly comparing the actions that generated them. 

Starting with the SU(2) pure gauge theory, where the gauge coupling is the only 
relevant parameter, the procedure is as follows. 

1. Generate an ensemble of gauge configurations with an action S{(3) on a lattice. 

2. Block these n times to produce an ensemble of configurations on a {L/s^Y lattice, 
and measure the expectation values of various observables on them. 

3. Generate a new ensemble of gauge configurations with an action S{I3') on a (L/s)^ 
lattice, for a range of values of f3'. 

4. Block each of these n — 1 times to produce an ensemble of configurations on a {L/ s'^Y 
lattice, and measure the same observables for each (5' . 

5. Interpolate in such that each observable after taking n steps on the larger lattice 
agrees with the same observable after taking (n — 1) steps on the smaller lattice. 

6. Repeat for different n, e.g. for s = 2, L = 32, three values can be used: n = 2, 3, 4. 

We have now identified, for each n, a pair of bare gauge couplings /?'), with lattice 
correlation lengths that differ by a factor s, |' = |/s. In the limit n ^ oo, the quantity 



is the step scaling function for the bare gauge coupling. This is the analog of the 
Schrodinger Functional step scaling function for the renormalised coupling, s), and 
in the UV hmit where g"^ ^ — 2A^/^, there is a simple relation between the two: 



(2.1) 



Sb{/3;s) _ a{u,s) 
/3 ~ u 



- 1. 



(2.2) 
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Figure 1: The two-lattice matching procedure described in Sec. |2} Horizontal arrows 
represent RG blocking steps, and vertical arrows indicate matched lattices. In this case 
the matching is done after 2(1) and 3(2) blocking steps on the 16'^(8'^) lattices. 



We use the following s = 2 RG blocking transform, labelled ORIG following the 
naming convention used in Ref. 



2n+2pL,v 



(2.3) 



where U is the original gauge field on a lattice, V is the blocked gauge field on a 
{Ljl)* lattice, Proj is the projection back into the group of SU(2) matrices, and a is a 
free parameter, which can be varied to optimise the transformation. Changing a changes 
the location of the FP, and how quickly the RT is approached in a given number of steps. 
Ideally it should be chosen such that 

• All operators predict the same matching coupling between {n,n — 1) pairs for a 
given blocking step n. (Deviations are a measure of the systematic error from not 
being at exactly the same point along the RT.) 

• Consecutive blocking steps predict the same matching coupling, i.e. the coupling 
for which (n, n — 1) pairs agree should be the same for all n. (Deviations show that 
the RT is still being approached in the irrelevant directions^ 



We also use two different s = 2 RG blocking transforms constructed using hypercube 
(HYP) smeared links [H], labelled HYP and HYP2 HYP blocking uses a product of 



^If there existed other relevant couphngs that had not been precisely tuned to their critical values, 
these deviations could also indicate an unwanted flow in these relevant directions 



3 



HYP smeared links with smearing parameters (a, 0.6, 0.3), 

= W[U]2n,,W[U]2n+,,, (2.4) 

where U is the original gauge field and Vr[f/]2;,^ is the HYP smeared link at (x, /i), defined 
in App. |X} 

HYP2 blocking is the same except that HYP smearing is applied twice to each link 
with smearing parameters (a, 0.3, 0.3). 

KT' = W'[UU,,W'[UU^,,, (2.5) 

Ideally our results would be independent of the choice of blocking, so the use of three 
different blocking transforms allows us to check the systematic errors of the procedure 
and the distance from the RT. The HYP and HYP2 blocking transforms have also been 
empirically found to work better than the ORIG transform at strong coupling [19i| . 



3 Pure Gauge Results 

As an initial test of the method, matching in /3 between 32(16) and 16(8) lattices was 
performed using ~ 2000 SU(2) pure gauge configurations for each /3. We matched in the 
plaquette, the three six-link loops, and three 8-link loops, shown in Fig. |2j See App. [B] 
for details of these observables. 



(a) P{ti,u) 




Figure 2: The plaquette, 6-link and 8-link gauge observables used. Explicit definitions 
are given in App. [B] 

An example of the matching of the plaquette is shown in Fig. |3(a)[ The red, green 
and blue horizontal lines show the average plaquette on the 32^ lattice after 2, 3 and 4 
ORIG blocking steps respectively, at /3 = 3.0, a = 0.57. The interpolated red, green and 
blue points show the average plaquette on the 16^ lattice after 1, 2 and 3 blocking steps 
respectively, as a function of also at a = 0.57. The value of P' where the two red lines 
intersect gives the matching coupling for n = 2, similarly the green and blue lines give 
the matching coupling for n = 3 and 4. 
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This matching is repeated for each observable, and the spread of predicted matchings 
for each n gives a systematic error on the central matching value. The whole procedure 
is then repeated for various values of a, as shown in Fig. 3(b) , to find an optimal value 
of a where subsequent RG steps predict the same matching value. The intersection of 
the last two blocking steps gives a central value for Sfe(/3 = 3.0; s = 2), while the range of 
couplings for which any of the blocking steps intersect within errors gives the uncertainty 
on this central value. This was repeated for other values of /3, and also with the HYP and 
HYP2 blocking transforms. Fig. |4] shows the resulting step scaling of the bare coupling 
Shi determined using ORIG, HYP and HYP2 blocking on both 32(16) and 16(8) lattices, 
along with the 1-loop and 2-loop perturbative predictions. In the scaling region Sh agrees 
well with the perturbative prediction, and the agreement between the different blocking 
transforms and lattice sizes shows that finite size errors are small, at least for the pure 
gauge case. 



SU(2) Pure Gauge; |i=3.0; ci=0.57; 32(16) plaquette matching 



p-ji' vs a. 32(16) lattice matching, ^=3.0 
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(a) Plaquette Matching 
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(b) Qf-Optimisation 



Figure 3: An example of the matching of the plaquette in /3 for the pure gauge case 
using ORIG blocking on 32(16) lattices. This is repeated for each observable to give a 
systematic error for each matching, then a is varied such that all blocking steps predict 
the same matching. 
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Figure 4: Bare step scaling Sb for the pure gauge theory. Matching is performed on 32(16) 
and 16(8) lattices using ORIG, HYP and HYP2 blocking transforms. All give results 
consistent with each other and with perturbation theory. 

4 Anomalous Dimension Results 

Having confirmed that the two-lattice MCRG method works for the SU(2) pure gauge 
case, we now turn to the full Minimal Walking Technicolor theory of two Dirac fermions 
in the adjoint representation of SU(2). 

There are two couplings of interest in this theory, the gauge coupling and the mass. 
At an IRFP the gauge coupling is expected to be irrelevant, leaving the mass as the only 
relevant operator. So we should in principle be able to match in the mass at arbitrary 
couplings, as long as we have sufficient RG steps for the gauge coupling to flow to its FP 
value, and the mass is small enough that despite the resulting flow in the mass the system 
remains close to the IRFP. In practice we only have a small number of RG steps, and 
because the beta function is small the coupling only flows slowly towards its FP value, 
so we set (3' = (3 assuming the flow in the coupling due to one RG step is negligible. The 
case (3' (3 will be considered later. We match observables as for the pure gauge case, 
but instead of matching in (3, we fix (3' = (3, and match pairs of bare masses (amo, a'Tn'o). 
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We use the HiRep [22] implementation of the Wilson plaquette gauge action with 
adjoint Wilson fermions and an RHMC algorithm with two pseudofermions. We generated 
~ 3000 configurations on 16^ and 8^ lattices, for a range of bare masses at each /3. This 
allows two matching steps, after 2(1) and 3(2) steps on the 16^(8'^) lattices respectively. 
See Fig. [5] for an example of this matching and subsequent a-optimisation. 

Because the bare mass is additively renormalised we convert the bare masses to PC AC 
masses, 

am = ampcAc ( 2 j ' «"^PCAc(a;o) = , — (4-6) 

where is the symmetric lattice derivative operator, and j a and /p are axial and 
pseudoscalar correlators respectively. We measure the PCAC mass, am, as a function of 
bare mass, amo, for each /3 on the 16"^ lattices, as shown in Fig. |6| We then use this to 
convert the bare masses on both 8^ and 16^ lattices to PCAC masses, as the measured 
PCAC masses on the 8^^ lattices suffer from finite volume effects, as described in App. [Cj 
Our previous result [23] for the anomalous mass dimension used PCAC masses measured 
on the 8^ lattices and hence contained a large finite volume effect, which has been removed 
in the present work. The anomalous mass dimension appears in the RG equation for the 
mass 

— — —r = -y^am = -(1 + 7)am. (4.7) 
dm 

At an IRFP the anomalous mass dimension is a constant, so the expression can be inte- 
grated to give 

a'm' ^,,1 , , 

= 2^+^ 4.8 

am 

for a pair of matching masses (am, a'm'), from which a value for 7 can be extracted. We 
used four values of (3, (3 = 2.15, 2.25, 2.35, 2.50, and the matching PCAC mass pairs using 
the HYP blocking tranform are shown in Fig. [7] We also repeated the matching using 
ORIG and HYP2 blocking. A comparison of the three blocking methods is shown for each 
(3 value in Fig. [9} and the agreement between them is very good. 

Different (3 values seem to predict consistent values for the anomalous mass dimension, 
as shown in Fig. [7} which uses all the beta values and masses in the range 0.02 < am < 
0.16. Masses above this range are excluded because the data are no longer consistent 
with a linear fit, indicating that we are beyond the small-mass linear scaling region. The 
masses below am = 0.02 are excluded because they are likely to contain a significant 
unwanted contribution from the running of the coupling, as described in App. |Dj 

The x^/d.o.f. for a linear fit of the form in Eq. 4.8 is shown in Fig. pi which favours a 



vanishing anomalous mass dimension. A combined fit to all /3 and m gives 7 = —0.03(13), 
with 7 > 0.13 ruled out at a 95% confidence level. Significant systematic errors exist that 
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have not been accounted for in this result, in particular the assumption that we can set 
/3' = /3. These systematic errors will be discussed in Sec. |6] 
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SU(2) Adjoinl; ^=2.50; a=0.65; am(,=-1 .05; 16(8) plaquette matching 



ifl' vs a, 1 6(8) lattice matctiing, p=2.50,am(,=-1 .05 




(a) Plaquette Matching 



(b) a-Optimisation 



Figure 5: An example of the matching of the plaquette in the bare mass amo. This is 
repeated for each observable to give a systematic error for each matching, then a is varied 
such that all blocking steps predict the same matching. 
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Figure 6: PCAC mass, am, as a function of the bare mass, arriQ, on 16^ lattices for 
(3 = 2.15,2.25,2.35,2.50 
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Figure 7: HYP Matching in mass using all /3 values in the mass range 0.02 < am < 0.16. 
Consistent with a vanishing anomalous mass dimension, 7 = 1 is strongly disfavoured. 
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X"/d.o.f. for HYP matching in range 0.02<am<0.16 
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Figure 8: x^/d.o.f. of 7 with HYP Matching in mass using all /3 values in the mass range 
0.02 <am< 0.16. Combined best fit with 22 d.o.f. gives 7 = -0.03(13). 
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Figure 9: Comparison of matching mass pairs using different blocking transforms. 
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5 Discrete Beta Function Results 



It would also be interesting to directly measure the step scaling of the bare coupling 
Sb{(3, s) for the full theory, as was done for the pure gauge case. A fixed point would be 
indicated by a change of sign in this quantity as the bare coupling is varied from weak to 
strong coupling. The difficulty is that the mass is a relevant operator, while the coupling 
is expected to be at best nearly marginal, so in order for the MCRG to pick out the 
behaviour of the coupling the mass would have to be tuned to zero. Furthermore, even 
if the mass is tuned sufficiently close to zero that we are initially following the evolution 
of the gauge coupling (assuming it is the least irrelevant remaining operator), we can no 
longer take n — )■ oo limit, because the coupling will flow to its FP value, and the flow of 
the mass will eventually dominate. 

So the method is to tune the mass close to zero initially, then take a few RG steps, 
where the flow is hopefully following the gauge coupling, and extract the running of the 
coupling from this, before the flow in the mass becomes signiflcant. 

We measure the PCAC mass on 16^ lattices for a range of /3 values at small masses, 
and for each /3 extract the critical bare mass using a linear interpolation in the PCAC 



masses, as shown in Fig. 10 



Critical bare mass measurements 
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Figure 10: Measurement of critical mo values from 16^ lattices. 



We then flt an interpolating function in (3 to these critical masses, as shown in Fig. 11 



which gives us the critical bare mass for any /3 in the range 2.15 < /3 < 2.75. The 
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systematic error in the PCAC mass due to this interpolation is small, < 0.001, as described 
in App. M 



Critical bare rrass interpolation 
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Figure 11: Interpolation of critical arriQ values from 16^ lattices, to be used to determine 
critical amg values for the massless 8"^ runs. 



The measured PCAC masses for the 16^ critical runs are shown in Fig. 12, along with 
some runs at small masses to allow us to check for any systematic mass-dependence. 

We generated ~ 3000 configurations on 16^ and 8"^ lattices for a range of (3 values, each 
run at the critical bare mass. The matching procedure in /3 is then essentially the same 
as for the pure gauge case. An example of the matching in the plaquette, and subsequent 



a-optimisation is shown in Fig. 13 



The resulting measurement of Sb{l3) is shown in Fig. 14 This includes both the 



massless and small mass 16^ runs - within errors Sb shows no mass dependence for these 
small masses. The ORIG matching values of Sb are clearly positive throughout, the 
HYP values are lower, and the HYP2 values are consistent with zero within error bars. 
There is no clear cross-over from positive to negative values of Sb for any of the blocking 
transforms, so while the data are consistent with a fixed point, they are not sufficiently 
precise to distinguish slow running from a fixed point. This level of precision is nonetheless 
similar to that found in the Schrodinger Functional studies, albeit with less understood 
systematic errors. 
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Critical mass measurements 
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Figure 12: Measured PCAC mass of critical 16^ runs used for massless matching 
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Figure 13: An example of the matching of the plaquette in /3' for the massless case using 
ORIG blocking. This is repeated for each observable to give a systematic error for each 
matching, then a is varied such that all blocking steps predict the same matching. 
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Figure 14: Massless matching in /3 on 16(8) lattices. This includes both the massless and 
some small mass runs at each /3 - within errors Sb shows no mass dependence for these 
small masses. 
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6 Systematic Errors 



6.1 Matching Observables 



In principle, for a given (/3, am) there should be a unique matching set of couplings 
(/3',a'm'), for which all blocked observables agree after n and [n — 1) blocking steps 
respectively. In this work we have set P = P' and we were able to find a'm' such that the 
blocked observables matched. 

In practice however, all of our observables are small Wilson loops, and as such are 
strongly correlated and have a very similar dependence on and a'm'. This means 
that we can in fact find a "matching" a'm' for a range of values of /?', which, given that 
we do not know the correct value of /?' to use, significantly increases the error on our 
determination of 7. As an example the matching mass pairs for /3 = 2.25 and various 



values of l3' are shown in Fig. 15 



Matching PCAC mass pairs - beta=2.25 
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Figure 15: Mass matching pairs at /3 = 2.25 for a range of /3'. Varying /3' can lead to very 
different results for 7. The dotted horizontal and vertical lines show mL, m'L' = 1. 



In Sec. [5} while we find that Sf, = /3 — /3' is compatible with zero, corresponding to 
setting /3 = f3', the error bars are relatively large, enclosing the region —0.08 ^ /9 — /9' < 
0.16. From Fig. 15 we see that for /3 = 2.25 this region is approximately bounded by 
f3' = 2.15 and /?' = 2.35, and encloses a large range of values for the anomalous mass 
dimension, —0.6 ^ 7 ^ 0.6. This range is representative of the errors in the anomalous 
mass dimension due to the uncertainty in the correct value of /?', and is the dominant 
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source of systematic uncertainty in our results. 

To resolve this issue, more observables which are "orthogonal" to the Wilson loops, 
for example meson correlators or other fermionic observables such as the PCAC mass, 
would need to be included in the matching determination to determine which of these 
blocked configurations are actually matched and find a unique matching set of couplings. 
We were unable to consider these observables with our current lattice volumes due to the 
very small size (2^ and 4^) of the blocked configurations on which the observables are 
measured. 



6.2 Finite— volume Effects 



The MCRG techniques that have been described previously have achieved much success 
when applied to QCD-like theories. However, previous studies both analytical and numer- 
ical of the minimal walking model studied in this paper suggest that it is likely conformal 
or very nearly conformal in the infrared. Thus the calculations of 7 require us to compute 
the flow in the mass parameter for a conformal field theory (CFT) in the presence of a 
small mass deformation. 

Of course to be a true CFT the theory must be considered in infinite volume but 
unfortunately when we do simulations we have to live with lattices of finite size and this 
makes the analysis a little more subtle [21]. In principle both the mass parameter and 
the lattice size then determine any correlation length. If we want to extract the correct 
physics for the infinite volume CFT it is important to make sure that the correlation 
lengths we are measuring and matching are not being strongly influenced by the finite 
box size. So we should make sure that the physical lattice size is much larger than the 
correlation length or equivalently m ^ 1/L where L is the lattice length, to obtain the 



separation of scales shown in Fig 16 




IRFP 



1/L m 1/a' 1/a UVFP 

Figure 16: Separation of scales in a normal lattice simulation. The lattice spacing a 
provides a UV-cutoff 1/a, the lattice size L provides an IR-cutoff 1/L, and the mass of 
the lightest propagating state m lies somewhere between these two scales. After an RG 
blocking step the lattice spacing changes to a' = 2a, and the modes between 1/a' and 1/a 
are integrated out. 
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Thus we are forced to consider masses that are sufficiently large to satisfy this con- 
straint. Unfortunately, if we use too large a mass we will move the system a long way 
from any IRFP and the simple MCRG techniques we are using will not apply. The only 
obvious way to reconcile these two things is to use large enough boxes that we can keep 
L ^ 1/m while making m small enough to keep us close enough to the CFT fixed point. 

For our L=16 lattices to achieve for example mL > 2 we would need am > 0.12, which 
may be simply too big to keep close the system close to the fixed point. 



7 Stability Matrix Method 

The two-lattice matching technique used in this work was first used to investigate quenched 
QCD EZ], and more recently QCD with many flavours of fermions [TSl [T2], and 

allows a determination of the most relevant coupling in a system. The original MCRG 
method [20] in principle allows the extraction of all critical exponents of a system, both 
relevant and irrelevant. 

7.1 Method 

Consider a hamiltonian that can be written as a sum of couplings Ki and observables Si, 

H = J2k^S,, (7.9) 

i 

and an RG transform Rs of scale s such that 

i/^'^+i) = = Kt^'^St+'\ (7.10) 

i 

where S^"''^^^ is the same observable as S^"'^ only measured on the lattice blocked n + 1 
rather than n times. The fixed point of the RG transform is defined by the condition 

H* = RsH* = J2 K*S*, (7.11) 

i 

and near this point the flow in the couplings can be expanded linearly to give 

^(n+l) _ ^ _ (7_^2) 



where 



(7.13) 

H* 
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The chain rule gives 

From which T^j can be constructed using the identities 

1^ = {St^S^-'^) - ^ 4) (7.15) 



= (^H^W) - {St^){sf) ^ (7.16) 

The eigenvalues of T*j give the critical exponents of the system [28] , e.g. z/ = In s/ In \h, 
where Xh is the largest eigenvalue, and so ?/m = l/j^ = In A/^/ In s. From a single simulation 
close to the critical point, correlation functions of blocked observables are measured to 
construct the matrix Tij, from which and other exponents can be determined. 

If we are sufficiently close to a fixed point, then the largest eigenvalue of T should stay 
constant as the number of blocking steps is varied, and also as the number of observables 
used to construct T is varied. This method requires a larger lattice and higher statistics 
than the two-lattice method, but potentially allows more information to be extracted, in 
addition to being a useful consistency check of the two-lattice method results. 



7.2 Pure Gauge Results 

We have seven blocked observables and four blocking steps on the 32^ lattices. This means 
we can vary the number of observables, and hence the size of the stability matrix T, from 
1 to 7. We can calculate T after 1/2, 2/3 and 3/4 blocking steps, for any choice of our 
blocking parameter a. Unlike in the two-lattice method there is no cancellation of finite 
size effects, so these are likely to be large. 

We applied this method to our pure gauge configurations, where we expect a marginal 
largest eigenvalue of 1. Using more than four observables (i.e. including 8- link loops) 
tends to give a complex largest eigenvalue of T, so we only use 1 to 4 observables to 
construct T. Fig. 17 shows the largest eigenvalue after 1/2 (red), 2/3 (green) and 3/4 
(blue) blocking steps, as a function of alpha using ORIG blocking, on a 32^ lattice at 
(3 = 3.0: 

The finite size effects are presumably very large, but in general this looks sensible - 
the variation with alpha is reduced as the number of blocking steps are increased, and 
for small a the eigenvalues are independent of the number of blocking steps within the 
spread of eigenvalues from varying the number of observables, and are consistent with a 
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Figure 17: The largest eigenvalue of the stability matrix T as a function of a, after 1/2 
(red), 2/3 (green) and 3/4 (blue) ORIG blocking steps on a 32'* pure gauge lattice at 
13 = 3.0. At each blocking step T is constructed using 1,2,3 and 4 observables, so the 
spread of eigenvalues of the same colour (same blocking level) gives some idea of the 
systematic variation from the number of observables included. 



marginal eigenvalue of 1 corresponding to the expected logarithmic flow of the coupling 
in quenched Yang-Mills. 



7.3 MWT Results 

For Minimal Walking Technicolor we have 16^ lattices so we can construct T after 1/2 



and 2/3 blocking steps, although again the finite size effects are likely to be large. Fig. 18 
shows the largest eigenvalue of T after 1/2 (red) and 2/3 (green) ORIG blocking steps 
for (3 = 2.25, am ~ 0.2. Between 1 and 7 observables are used to construct T, and the 
spread of eigenvalues at a given blocking level is small, showing little dependence on the 
number of observables used. On the other hand there is a large difference between the 
two blocking steps, which suggests that we are not close to a fixed point. 

This is representative of the situtation for all of our runs - the picture is qualita- 
tively the same for virtually all of our values of (3 and m, and for all three RG blocking 
transforms. 

These results suggest that while the pure gauge simulations are sufficiently close to the 
renormalised trajectory, the MWT simulations are not close enough, and more blocking 
steps, and hence larger lattices, are required. 
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Figure 18: The largest eigenvalue of the stability matrix T as a function of a, after 1/2 
(red) and 2/3 (green) ORIG blocking steps on a 16^ lattice at /3 = 2.25, am ^ 0.2. At 
each blocking step T is constructed using between 1 and 7 observables, so the spread gives 
some idea of the systematic variation from the number of observables included. 

8 Conclusion 

We find a vanishing anomalous dimension and slow running of the coupling consistent 
with an IRFP, however our results suffer from considerable systematic errors. The most 
significant of these is the uncertainty in the location of the fixed point in the coupling, 
which in turn produces a large uncertainty in the value of the anomalous mass dimension. 
This is essentially the same issue that the Schrodinger Functional studies of this theory 
have suffered from. 

Assuming negligible running in the coupling we find 7 = —0.03(13). Using our errors in 
the measurement of the discrete beta function to suggest a window for the possible flow in 
beta produces a large uncertainty in the associated anomalous dimension —0.6 ^ 7 < 0.6. 
Previous studies |S1 El El E] have provided evidence of a backward flow in the coupling 
at these values of (3 consistent with an IRFP. If this is the case our results indicate that 
the estimate of 7 = —0.03(13) assuming no flow in (3 would actually constitute a lower 
bound on the true anomalous dimension. A small value of 7 is also likely to result from 
finite volume effects; if the beta function is small but non-zero then many blocking steps 
and hence large starting volumes will be needed to move the system away from the usual 
perturbative UV fixed point. 

While the MCRG method is potentially a promising technique for studying possible 
conformal systems our simulations indicate that it is currently limited by several sources 
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of systematic error. Perhaps the single largest factor contributing to this error is the 
relatively small lattices that have been used in this study; this places constraints on how 
small a quark mass can be used and limits the number of blocking steps that can be 
taken. These restrictions mean that in all likelihood the RG trajectory does not closely 
approach the putative IRFP making extraction of an exponent like 7 problematic. 

Adding more matching observables, in particular fermionic ones such as meson corre- 
lators, should give a broader set of observables which will enable us to determine unam- 
biguously matched actions. Using 32^ and 16"^ lattices matching could be done after 1/2 
and 2/3 blocking steps - the smallest blocked lattice being of size 4^, so that for example 
the meson correlator at each timeslice could be used as an observable. 

Going to 32"^ lattices will also allow us to go to smaller masses and thus be closer to 
the IRFP, as well as allowing us to match our current set of blocked observables after 
three blocking steps instead of two, which is a much more stringent matching condition 
that should help to constrain /?'. With larger lattices we will also be able to check for 
finite volume effects using the analysis of Ref. [29]. In addition it may be possible to 
get a determination of the anomalous mass dimension using the stability matrix MCRG 
method. Other potential improvements include using an improved action to reduce 0{a) 
effects, and adding an adjoint plaquette term to the action to move the system away 
from a non-physical ultraviolet fixed point due to lattice artefacts, and allow us to go to 
stronger coupling [29] . 

However, even taking into account our currently large systematic errors, we find a 
small anomalous dimension which is clearly less than 1, and hence that, in its simplest 
form [30], the SU(2) gauge theory with two Dirac fermions in the adjoint representation 
is not a viable Walking Technicolor candidate. 
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A HYP Smearing 

A HYP smeared link is constructed in three stages [2T], each stage consisting of 

a modified APE blocking step [3l] projected back into the gauge group. In the first stage 
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a set of decorated links V are constructed from the original links U, 



( \ — ^ t 



Prom these a second set of decorated links V are constructed, 



(A.17) 



Vn,p.-v = Proj 



U-2 — ^ — ~ \ 



(A.18) 



and finally the HYP smeared hnks are constructed from this second set of dressed links, 



W[U] 



n,p 



Proj 



^1 - ai)Un,p. + Vn,u;pyn+u,p.;uV^ 



t 



B Observables 



Explicit expressions for the seven observables used in the matching: 



Ej=i T.u=o ^6(/^, i^, p = (/X + 1) mod 4) 



(/i + 2) mod 4) 
{ji + 3) mod 4) 
(/i + 1) mod 4, a 



(/i + 1) mod 4) 



E^=i Y.y=o ^8(/^, P = (/^ + 2) mod 4, a = (/i + 1) mod 4) 
Ej=i El^=o ^8(p, i^, p = (p + 3) mod 4, a = (p + 1) mod 4) 



where 



Lsil^,^^, p) 



E.^Tr 
E.^RTr 
E.J^Tr 



Ux,uUx-\-u,uU,_ 



ul 



'c+v+p,p^ x+v,p^ x,v 



(A. 19) 



(B.20) 



Ux,pUx-\-p,^uUx+fi+v,pU 

Ux,pUx+fi,uUx+fi+l/,pUx+p,+u+p,aUl_^_,y^p_^_a,p.Ul-^-l/-^-p^a^x+l/,p^x,u 

(B.21) 



C PCAC masses 

The PCAC mass measured on the 8^ lattices suffers from a finite- volume effect, the time 
extent of the lattice is not sufficient for the mass to reach a plateau. This can be seen in 
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Fig. 19, where the PCAC mass for the same values of {P,am) on lattices of size 8^ and 
16^ are compared. 

The lattice artefacts however seem to be small - the PCAC mass as a function of t 
agrees well between the two lattices. We use the PCAC mass measured at t = 8 on the 
16^ lattices throughout this paper to convert bare masses to PCAC masses. To estimate 
the systematic error due to this procedure, we can compare the PCAC mass at the same 
timeslice t = 4 on the two lattice sizes, as shown in Fig. 20 The difference is smaller than 
the statistical error on the PCAC mass for small masses, and increases approximately 
linearly with mass to ~ 0.005 at am = 0.20. To include this source of error, a systematic 
error of 0.025am was added linearly to the measured error on each matching mass. 



D Matching mass range 

We have assumed that the contribution from the irrelevant directions is negligible by the 
time we perform the matching. But in reality we match after only one or two steps, and 
there may still be a contribution coming from the coupling. We are forcing the matching 
coupling to the same value as the initial coupling. This causes a shift in the measured 
observable, which is then cancelled by a shift in the observed matching mass from the 
true value. 

Consider the situation where we measure some observable C»(")(/3, am) on the larger 
lattice after n blocking steps. The correct matching done in both the coupling and the 
mass would give a'm', /5' such that 

C("-^)(/3',aW) = 0^''\(3,am) 

Instead we fix = /3 and find a"m" such that 

Dropping the ^^~^') superscript on O and taylor expanding around the correct matching 
values gives 

a"m") = 0{P', a'm') + cp{(5 - (5') + c^{a"m" - a'm') 

where 

_ dO dO 

This gives the relation 



/3',a'm' 



II II II I Q Ql\ 

am = am [p — p ). 

Cm. 
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So the effect of not matching in P' will cause the measured a"m" matched to am = to 
be shifted away from zero. 

For this reason we excluded the matchings at am < 0.02 from the anomalous mass 
dimension fits, since they are likely to have significant contributions from the running of 
the coupling, although including them makes little difference to the fits. 

In fact the size of this shift at am = could give us an indirect way to measure the 
running of the coupling. At am = 0, we know that a'm' = 0, which gives the relation 

{13-13') = -^a"m". 

The difficulty in determining Cm and c^, as well as the relatively large uncertainty on the 
measurement of a"m", means that although in principle this would be a way to measure 
(/3 — /?'), in practice it is not feasible with our data. 



E Massless interpolation 

The difference between the measured critical bare masses, and the values predicted by 
the interpolation function, are very small as shown in Fig. 21 This implies a systematic 
error < 0.001 on the PCAC masses due to the interpolating function. 
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PCAC mass vs t; p=2.25; mo=-1 . 1 
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gure 19: PCAC mass as a function of time on 8^ and 16^ lattices. 
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Comparison of PCAC mass at t=4 between L=16 and L=8; |3=2.25 
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Figure 20: PCAC mass at t=4 on 8^ and 16^ lattices. 
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Difference between measured critical bare masses and interpolated values 
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Figure 21: To estimate the systematic error, here is the difference between measured 
critical bare masses, and the values given by a 3-parameter interpolation function fitted 
to only three of the measured values. Also shown is the difference between measured 
values and a 4-parameter interpolation function fit to all the data. The difference in 
PCAC mass is (1.5 — 3.5) x this bare mass difference, which is still < 0.001, and smaller 
than the statistical uncertainty on the measured PCAC mass. 
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